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We solve exactly the non-Markovian dynamics of a cavity mode in the presence of a feedback 
loop based on homodyne measurements, in the case of a non-zero feedback delay time. With an 
appropriate choice of the feedback parameters, this scheme is able to significantly increase the 
decoherence time of the cavity mode, even for delay times not much smaller than the decoherence 
time itself. 
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^. ■ Although feedback schemes have been used for long times to control noise, a general theory of feedback for quantum 

p^j \ systems has been developed only some years ago by Wiseman and Milburn [&p|. Interesting possibilities are opened 
by the ability to control systems at the quantum level using appropriate feedback loops and some of them have been 

t— I ' showed in a series of papers ||-|J. Ref. [Q has shown that an electro-optical feedback loop based on homodyne 
measurements of a cavity mode provides an affordable way to realize a squeezed bath for the mode. As a consequence, 
homodyne- mediated feedback can be used to get squeezing [p| , and, in the case of optical cavities with an oscillating 

L^ \ mirror, it can be used to significantly cool the mirror. This fact can be extremely useful for the interferometric 
detection of gravitational waves || ■ The application of a feedback loop realizes an effective "reservoir engineering" 
|l0f and therefore it can be useful also for decoherence control, which is a rapidly expanding field since decoherence 
is the main limiting factor for quantum information processing |0|. Refs. M ffl|| have already shown that the 

0^ ' decoherence induced by photon leakage in electromagnetic cavities can be significantly suppressed with appropriate 
feedback loops, using the homodyne photocurrent in |||| and direct photodetection and atomic injection in |j],|j. 

However, all the relevant applications considered up to now always assume the zero feedback delay time limit r — > 0, 
which is much easier to handle because the problem becomes Markovian and the effect of feedback can be expressed 
<—j ■ in terms of an effective master equation 0-0] . The presence of a non-zero delay has been considered briefly only in 
P| , where the spectrum of a homodyne measurement has been evaluated for a simple case. The Markovian treatment 
is justified whenever the feedback delay time is much smaller than the typical cavity timescale. If one considers 
squeezing or some other stationary state phenomenon, the feedback delay time r has to be compared with the cavity 
relaxation time 7 -1 and for sufficiently good cavities the Markovian condition 7T <C 1 is usually satisfied. However, 
if one considers the feedback scheme for decoherence control, the delay r has to be negligible with respect to the 
decoherence time tdec — (7"-) _1 j which can be much shorter than the damping time when the cavity mean photon 
C^ ■ number n is large. In these cases, the unavoidable non-zero feedback delay time may have important effects and it 
would be important to deal with the exact non-Markovian problem with t/0. There is in fact a renewed interest 
in non-Markovian effects, which can play an important role when considering quantum optics in high-Q cavities and 
in photonic bandgap materials. For this reason, non-Markovian trajectory theories have been recently developed in 
Refs. ff-|l4 
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uantum theory of feedback has been developed by Wiseman and Milburn in pl,pf using quantum trajectory 
15[ , and only later Wiseman showed an equivalent derivation based on the input-output theory fllq-Jlq] in 
Ref. U . However Ref. proved the equivalence between the two approaches in the perfect detection 77 = 1 case only. 
In this paper we shall see how to extend the quantum Langevin approach of the input-output theory to the non-unit 
efficiency case and we shall see that this theoretical framework is best suited to deal with the non-Markovian case of 
non-zero feedback delay time. We shall consider the non-Markovian effects by completely solving the dynamics of a 
cavity mode in the presence of a homodyne-mediated electro-optical feedback loop, which has been already considered 
(in the zero delay limit only) in Ref. M . 

The paper is organized as follows. In Section II we shall reconsider the quantum theory of feedback in the case 
of homodyne measurements, adopting the input-output theory of Gardiner and Collett p6|dl8| and we shall see how 
to introduce the non-unit detection efficiency in this framework. In Section III we shall completely solve the non- 
Markovian dynamics in the presence of a non-zero feedback delay by considering the time evolution of the probability 



distribution of the measured field quadrature and of the characteristic function. We shall consider in particular the 
possibility of inhibiting the decoherence of a Schrodinger cat state initially generated in the cavity and we shall see 
that the significant decoherence suppression which can be obtained, for appropriately chosen feedback parameters, in 
the zero delay case (see Ref. ||) is recovered even in the presence of not negligible feedback delay times. 

II. HOMODYNE-MEDIATED QUANTUM FEEDBACK THEORY WITHIN THE INPUT-OUTPUT 

FORMALISM 

We shall consider an optical cavity, with annihilation operator a, subject to the homodyne measurement of the field 
quadrature 

X lp = hae- iv + aU^). (2.1) 

We shall consider the possibility of applying a feedback loop to this cavity mode, by feeding back part of the output 
homodyne photocurrent to control in some way the a mode dynamics. 



using the input- 



First of all it is convenient to reformulate Wiseman and Milburn quantum theory of feedback p],p 
output theory developed by Gardiner and Collett |16|Jl 7| ] . The input-output formalism is essentially an Heisenberg 
approach for the whole system (cavity and vacuum bath), in which the environment dynamics is described by the 
white-noise input operator dB{t) satisfing the Ito rules |H] pj| : 

dB{tf = dB^(t) 2 = 0, dB^(t)dB(t) = 0, 

(2.2) 
dB(t)dB*(t) =dt, 

and the following commutation' s relations: 

[dB(t),dB(t')] =0, 

(2.3) 
[dB(t),dB^(t')] =6{t-t')dtdt'. 

In the absence of any feedback loop, the evolution of a generic cavity mode operator 0(t) in the interaction picture 
is described by the quantum Langevin equation [|16|Jl7| : 



7 I 
V7 [0(t), o+(t)] dB(t) + ^7 [6(t), a(t)] dBHt), (2.4) 



d6(t) = -(2a){t)6{t)a{t) - a) '(t)a(t)6(t) - 6{t)a\t)a{t) ) dt 



where 7 is the cavity damping rate. 



Eq. (2.4) can be solved explicitcly in terms of the evolution operator U(t,to) 

d(t) = U\t,0)OU(t,0), (2.5) 

which in the absence of feedback takes the following form Eq| 

U(t,to) =«P (V7 / dB\t')a - V7 / dBfflat), (2.6) 



to 



where a and a^ are in the Schrodinger representation and exp denotes the time-ordered exponential. The evolution 
operator U(t,t ) describes also the evolution in the Schrodinger representation, |18| , |l1| , 

|*(*))r = !7(* > 0){|V>4)®|0)}, (2.7) 

where the vector \^(t))T obeys the following stochastic equation 

d|*(t)) r = {V^a dBHt) - V7a f dB(t) - | a t adi}|*(t)) T . (2.8) 



Using the commutation' s rules (2.3), it is easy to prove that the Heisenberg evolution (2.5) satisfies the usual 



requirement that the input noise dB(t) has to commute with every cavity operator evaluated at preceding times t' < t 

[dB{t),O(t , )] = [dB^{t),d(t')]=0 for t>t'. (2.9) 

Equation fl2.4| ) can be used to get the time evolution of a generic matrix element of 0(t) between two state vectors 
of the whole system of the form \4>a) ® |0) and \iPb) <8> |0), in which the environment is left in the vacuum state, 

d(0(t)) AB = 1 { 2(aHt)6(t)a(t)) AS - (o t (t)o(t)6(t)) AB - (0(t)a!(t)a(t)) AB } dt, (2.10) 



where 



(6(t)) AB = {<o| ® (i> A \}d(t) {\ii> B ) ® |o)}. (2.ii) 



Let us now introduce the feedback loop associated to the homodyne measurement of the quadrature X v . Differently 
from Ref. |3J, we assume the possibility of a non-unit homodyne detection efficiency rj < 1. The application of a 
feedback loop is equ ivalent to add a feedback Hamiltonian Hfb(t) 0J3J§) so that the correction to the Heisenberg 



evolution of Eq. (2.4) takes the form 

d d fb (t) ^^ dY ^- T) [F { tld(t)] 



V 



(2.12) 



where F(t) is the observable of the cavity mode through which the feedback acts on the system and Y v (i) is the 
output field operator associated to the homodyne measurement. In the quantum trajectory approach of [0,p|, the 
fed-back homodyne photocurrent Y v (t) is a classical quantity, but in the quantum Langevin approach it must be an 
operator with its quantum fluctuations. However, one can adopt the general theory of homodyne measurements of 
Ref. p0[ and write the photocurrent operator in an analogous way 

dY lf> {t)^2^r,X v {t)dt + v fT 1 d~ v (t). (2.13) 

where dS v (t) describes the "noisy" part of the output photocurrent operator. The only delicate point in the derivation 
of the quantum theory of feedback of pLpj in the imperfect detection case using the input-output theory is just 
the exact det ermin ation of this noisy operator daJtl Since Y v (t) is an output operator, it is quite natural to 
consider Eq. ([2.13|) as an input-output relation p6 17 19], so that the nois y term would simply be the input field 



dS ¥ ,(t) = dB{t)e~ lip + dB^{t)e %v . However this interpretation of Eq. ( ]2.13| ) is correct only in the perfect detection 
case r\ = 1, because only in this case the output of the detection apparatus coincides with the cavity output and 
the quantum fluctuations of the vacuum bath are transferred unaltered by the detector. In the presence of imperfect 
detection, the output photocurrent may be non trivially related with the input noise dB(t) and in general one has to 
describe the noisy operator dEL(i) in terms of a new noise dBf(t), which we shall call "feedback" noise. Therefore 
one has to write 

dE ip (t) = dB f (t)e- ilp + dB\{t)e lip , (2.14) 



where the feedback noise dBf(t) satisfies the same properties (2.2) and ( |2.3|) of the input noise; moreover this feedback 
noise is correlated with the input noise dB{t) and this correlation is determined just by the detection efficiency 77, 
since one has 

[dB{t),dB f (t')] =0, 

(2.15) 
[dB{t),dB){t')} = ^/TjS{t - t')dtdt' . 

It is immediate to see that in the perfect detection case 77 = 1, one can identify the feedback noise with the input 
noise dBf(t) = dB(t), while in the opposite case r\ = the two noises are uncorrelated, as it can be easily expected 
since in this case the fed-bac k noi se has nothing to do with the vacuum input noise. 



In the feedback correction ( 2.12 ) of the Heisenberg evolution, r is the delay time associated to the feedback loop and 
since it is a non- negative quantity, it ensures that the output operator Y v (t — r) commutes with all system operators 
evaluated at time t. In particular Y v (t — r) commutes with F{t) and so there is no ambiguity in the definition of 



dOfl(t). As it has been stressed in Ref. 0,§,||, one must be careful in using (2.12); the feedback process is physically 



added to the evolution of the system of interest, so its stochastic differential contribution has to be introduced as limit 



of a real process. This implies that (2.12) has to be considered in the Stratonovich sense. Therefore it is convenient 



to rewrite it in the Ito form and then add it to Eq. ( |2.4| ), so that the resulting equation for 0(t) becomes 
dO(t) = 2( 2at(t)d(t)a(t) - at(t)a(t)6(t) - 6(t)<J (t)a(t) ) dt 
-^[6(t),aHt)]dB(t) + ^Y[d(t),a(t)]d&(t) 

+i V7 ( ^{t - t) e*<> dt + d -^f±e i * ) [F(t), 6(t)] Q(t - r) (2.16) 

+i V7 [F(t),6(t)] ( Sya(t - r) e-** dt + d -El^ll e ^ ) @{t _ r ) 

-gjQ(t-T)[F(t),[F(t),6(t)]]dt. 

This is the quantum Ito stochastic equation describing the time evolution in the presence of feedback and non-unit 
detection efficiency and it coincides with the quantum Ito equation derived in Ref. p| in the case of perfect homodyne 
detection r\ = 1. We have also explicitely inserted the step function 0(i — r) with respect to Ref. [J3| to stress the 
impossibility for the feedback to act on the system before the delay time r has elapsed since the initial condition. 
This means that the evolution of 0(t) for < t < r coincides with that in absence of feedback, described by Eq. (^J). 
Moreover we observe that, since the equation for dO(t) contains only stochastic ter ms evaluated for times t' < t 
(precisely t and t — r), it is possible to con clude that the commutation relations (|2.9|) are valid also in the presence 



of feedback and, more in general, that Eq. (2.16) preserves the canonical commutation rules for a(t) and a'(t). 



A. Zero delay time limit 

Up to now, the explicit applications of the quantum theory of feedback of Wiseman and Milburn have considered 
the zero delay time case r = only, when one has a tractable Markovian equation. Whenever one considers a non-zero 
delay the problem becomes non-Markovian and difficult to solve. 

The feedback master equation for homodyne-mediated feedback in the zero delay limit has been first derived in its 
general form using quantum trajectory theory [|15| in Ref. [Qg]. In the case of perfect homodyne detection r\ = 1, the 
same homodyne-mediated feedback equation has been rederived using input-output theory by Wiseman in [pi and, 
in its linear stochastic form, by Goetsch et al. in Ref. [p|. However, the connection between this linear stochastic 
approach and the input-output theory was not made explicit there. In reviewing the zero delay time case we shall 
clarify here the connections between the different approaches and we show in particula r tha t the linear stochastic 
Schrodinger equation app roach of Ref. || is equivalent to the input-output resul t of Eq. ( 2.16J ) (in the case r\ = 1) in 
the same way as Eq. (|2.8|) is equivalent to the quantum Langevin equation Eq. (|2.4|). 



The starting point of the analysis of Ref. |6j is the evolution equation for the state vector \^>(t))T of the whole 
system (cavity and vacuum bath). In the no feedback case, this equ ation is obviously equal to Eq. ( [2.8| ); the feedback 
loop is then introduced using the same Hamiltonian modification of ( [2.12 ) in the case 77 = 1. However, since in Ref. 



the zero delay time limit is considered from the beginning, one has to be careful with operator ordering, because in 
this circumstance one is not guaranteed that rfY^(t) commutes with the cavity mode operator F(t) at the same time. 
In Ref. 6l the question is solved by imposing that "the feedback acts later", i.e., that \^(t + dt))x is obtained from 
|\&(£))t by means of U(t + dt, t) (the evolution operator in the absence of feedback loop) first and by the Hamiltonian 
feedback correction later (see Eq. (2.11) of B). It is clear that in this way the equivalence with the input-output 
approach of Ref. ||] it is not so evident. However the equivalence can be proved by considering the following evolution 
operator 



U F (t,t ) =^xp { j'l [y^dB^{t')a - ^dB{t')a) - i^jFe~^dB(t') 
-i^jFe itp dB\t') - i^Fe^dt' - qFae^ v dt^\; 



(2.17) 



where a, a? and F are in Schrodinger representation. Using the Ito rules (2.2) to evaluate the different ial d O(t), it 
is possible to check that UF(t,to) is just the evolution operator determining the formal solution of Eq. ( 2.16 ) in the 
t = and 77 = 1 limit, according to the usual rule 

6{t) = ul{t,0) OU F (t,0) . (2.18) 



As it can be easily expected, the zero delay evolution operator Upit,^) reduces to the no-feedback one U(t,t ) of 
Eq. ( [2.6| ) when F = 0. This suggests that an equivalent Schrodinger representation could be obtained also in the 



case of feedback with zero delay time, starting from an equation analog ous to Eq. (2.7). In fact if we apply £/i?(£,0) 



to the initial state of the total system and we use again the Ito rules (2^), one gets the following linear stochastic 
Schrodinger equation, 

d\^{t)) T = l^jadB\t) -Vra 1 dB(t) - -a)adt 

-iyfy F dE v (t) - ijF ae- ltp dt - | F 2 dt\\${t)) T (2.19) 

coinciding with the equation obtained in Q. This shows that the approach of Ref. |J and that of Ref. ||] are 
respectively the Schrodinger and Heisenberg view of the same theory, with C/p(i,to) the unitary operator mediating 
the transition from one to the other. 



In Ref. fdl, by adopting an appropriate representation basis for the vacuum modes (see for example [191), Eq. (2. IS) 
was then reduced to a linear stochastic equation for the cavity mode only which was solved numerically. In Appendix 
A we shall reconsider this linear stochastic equation for the cavity mode and we shall see how it is possible to solve 
it analytically by adopting the integration method described in [ pl[ . 

Determining an evolution operator analogous to Up(t, to) for the t^O case is much more difficult and we shall not 
consider this strategy to study the non-zero delay problem. Instead wc shall adopt the input-output formalism which 



has yielded Eq. (2.16). To be more specific, we shall alw ays co nsider generic matrix elements as those of Eq. (2.11) 
whose evolution equation can be easily derived from Eq. (|2.16): 



i(6(t)) AB = I { 2<ot(«)(5(*)o(t)) AB - (aHt)a(t)6(t)) AB - (d(t)a^(t)a(t)) AB } 

+ l7 {(« t (i-r)[f(t),6(t)]) AB ^ + ([F(t),0(t)]a(t-r)} AB e-^}&(t-T) (2.20) 

-%{[F(t),[F(t),6(t)]]) AB .Q(t-r). 

III. FEEDBACK DYNAMICS IN THE PRESENCE OF A NON-ZERO DELAY 

The dynamics in the presence of a feedback loop with a non-zero delay time has never been completely solved 
because of its intrinsic non-Markovian nature. In this paper we shall analyze the effects of a non-zero feedback delay 
by considering a specific example for the "feedback operator" F(t) 

F{t) = gX e (t) = |( o(«) e-' w + <**(*) e w ), (3.1) 

where the consta nt a represents the gain of the feedback process and is an experimentally controllable phase. The 



particular choice (3.1) of F means that the feedback loop adds a driving term to the mode dynamics, which could be 



achieved, e. g. , by using an electro-optic device with vari able transmittivity driven by the homodyne photocurrent 



The homodyne-mediated feedback model with the choice (3.1) for F(t) has been completely solved in Ref. |6J] in the 



Markovian limit of zero delay time and therefore the comparison with the results of Ref. [p| will be very instructive. 
As it is shown in Ref. B , the main virtue of the homodyne-mediated feedback is its capability of slowing down the 
decoherence associated with cavity damping provided that the feedback parameters g and are appropriately chosen. 
Here we shall see that the decoherence inhibition caused by the feedback takes place also in the presence of a non-zero 
feedback delay time; in particular, decoherence is appreciably slowed down even for delay times not much smaller 
than the decoherence time itself. 

First of all we shall show the exact time evolution for the marginal probability distribution P{x ip , t) of the quadrature 
component X^it): this will result in a quite simple expression which can be easily analysed. Then we shall give the 
complete solution of the system dynamics in terms of the symmetrically ordered characteristic function. 



A. The marginal probability distribution 

A good, even if not complete, description of the state of the cavity mode is given by the marginal probability 
distribution P(x v ,i) of the measured quadrature component X v (t). We shall consider the following class of initial 
states for the whole system 

PT (0H^A„,»(/3|®|0)(0|, (3.2) 

oc,0 

i.e. a linear superposition of coherent states for the cavity mode and the vacuum state for the electromagnetic bath. 
We shall focus on the evaluation of the moments 

(X»(t)) Pa = {<0| ® (0\}x£(t) {\a) ® |0)} (3.3) 



(N is an integer). P(x v ,0) is a superposition of Gaussian functions and, thanks to the choice of Eq. ( |3.l| ) for the 
feedback operator, the evolution in the presence of feedback remains linear, so that P(x ip ,t) will maintain its initial 
Gaussian behaviour. This fact will be explicitely verified at the end of the section. 

We now proceed step by step: first we explicitely determine the evolution of the first and second order moment and 
then we derive a recursive relation between the moments clearly showing the Gaussian nature of the corresponding 
probability distribution. 



For N = 1 Eq. (3.3) becomes the matrix element of the measured quadrature X v (t), and Eqs. (2.20) and (3.1) 



yield the following differential equation for its evolution 

j t (X v (t)) AB = -l(X v (t)) AB + 19 sm(9 - <p)(X v (t - r)) AB Q(t - r), (3.4) 

which can be integrated for t > using Laplace transforms as it is discussed in Appendix B: 

(*v>(t)>AB = (^(0))abX(*) (3-5) 



with 



Xd) = E ^V^" e- ¥t - nT) (? (« - -))"- (3-6) 



where [x] indicates the integer part of the real number x. In the case of the initial condition (|3.2[), one has simply to 
consider \ipA) = \P) and \i/jb) = \ot) in Eq. ( |3.5| ), eve n th ough it is clear that this solution is valid for any choice of the 



initial state of the cavity. The function %(t) of Eq. (3.6) will often appear in the complete analytical solution of the 
problem described in the following and it is therefore useful to describe its behavior in the various limits. The solution 
(3.5) has the correct behavior both in the g — ► limit, where one has the simple exponential decay x(t) = e~ 7 */ 2 , 



and in the r — > limit, yielding 

X (t) = exp { - | (1 - 2gsm(8 - <p)) t), (3.7) 

which is the same solution which can be derived from the exact treatment of Ref . . It is interesting to consider the 
limit of small delay, jt <C 1, because this condition can be easily verified using good optical cavities and common 
electro-optical feedback loops. From the exact solution ( |3.6| ) one gets 

x{t) = {! _ g^(g-y) ( 2 _ lt{l _ 2gsin{e _ ^ 7r j e - { ^ 9 sMo- v m K (3 . 8) 

We have plotted the solution ( |3.6| ) in Fig. H, in which there is also a comparison with the no feedback case and with 
the Markovian feedback case of the zero delay time limit. This plot makes evident how the major part of the difference 
between the zero delay and the delayed cases comes from the retardation caused by the presence of <d(t — r): for 
simplicity, we shall refer in the following to this effect as th e "s tep effect" . 



For the determination of the second order moment (Eq. (3J3) with N = 2) it is convenient to study the correlation 
function C(t,t') = (X v (t)X v (i/)) 0a for ev ery positive values of t and t' . We first focus on the dependence on t' and 



differentiate C(t,t') by t' using Eq. ( 2.16 ), so to get the following differential equation: 



dC(t,t') 



7 



C(t, if) dt' + 75 sin(0 - <p)C(t, t' - r)9(t' - r) dt' 



(3.9) 



-^e^<[^(*),d£+ (/)]>,„ + J^dMO ~ vy*{[X v {t),dB\{t' ~ r)]),Mt' r). 



The corresponding initial condition can be obtained from Eq. (3.5) and is given by C(t,0) = (X v (0)X v (Q)) p a x(f)- 
Now Eq. (pi]) can be explicitely integrated once that the explicit t' dependence of the commutators between the 
noises dB^(if) and dB\{if — r) and the operator X v (t) is known. For the determination of these commutators, let 

us define first of all f(t, t') — [X v (t), B^ [t')\ , in which B^(t') is the sum of all the Ito increments of the input noise 
from the initial time up to t' . Now we can proceed in an analogous way as we have done to write Eq.(3.£) and we 
differentiate in t by keeping t' constant. Using the commutation rules of Eq. (2.3) we get 



§- t f(t, = -\f(t, t') +7.9 sin(0 - <p)f(t - r, t')Q(t - t) 



(3.10) 



y/1. 



-lif 



|e(i' - i) -gsm{9 - ip)Q(t - r)9(t' -t + r)\. 



which, except for the presence of last term which is a known function of time, is an equation similar to Eq . (3.4) 
and therefore can be solved using Laplace transforms and the initial condition f(0,t') = implied by Eq. (2.9). The 
commutation rules between X v (t) and the Ito increment dB'(if) can now be obtained by simply differentiating this 
solution with respect to t' and the final result is 



[X v {t),dB\t')} = F(t,t')dt' 



where 



f(t,t') = -^le-^{0(t - t') X (t - t') - gsm(d - <p) G(t - t' - r) X (t -if -r))=T{t- t'). 



(3.11) 



(3.12) 



The same procedure can be adopted to determine the commutator [X v (t), dBl(t')\ involving the feedback noise and 
one finally gets 



[X v (t),dBUt')] =F f {t,t')dt' 



where 



ff(t,t') 



V7, 



-lip 



{^j&(t-t')x(t-t') 



gahx(6-ip) 



(t-t'-T) X {t-t' -r))=T f {t-t'). 



(3.13) 



(3.14) 



dB^(t'), dBf(t r ), dBl(t') and a(t), d^(t') which we explicitely report here because they will be useful in the following 



2 l v ' * — ' 0, 

Proceeding as before it is possible to compute also the commutation rules between the noise operators dB(t'), 

[a(t),dB(t')] -=K{t,t')dt' [o(«),dBt(f)] = -n\{t,if)dif 

[o+(t),dBt(f)] = -TZ* {t,t')dt\ 
[a(t),dB}(t')] = -Tl{*(t,if)dt' 
[ot(«),dBj(f)] =-1lf*{t,t')dt', 



[a*(t),dB(if)] =n 1 (t,t')dt' 
[a(t),dB f (t')] =Kf(t,t')dt' 
[at(t),dB/(f)] =1l{{t,t')dt! 



(3.15) 



where 



{t _ t ,_ T) X (t-t<)-e-l^') 
gsm(6-ip) 



<R(t, if) = CQ- ge l{e+ip) {<d{t -t'-r) 

+G(t -if - 2r) [x(t -t'-r)- e-i^-*'-^} } := K(t - t') 



(3.16) 



fti(t,0 = -2e- iv F*{t,t') - e- 2lv K(t,t') = K^t-t'). 



llf(t,t')=iJ^ge i « > +rt{9(t-t'-T) 



gsin(6- tp) 



(3.17) 
(3.18) 



+&{t -t'- 2r) [ X (t -t'-r)- e -7(*-t'-r)/2] \= n f(j._ ^ 



ft{(M') = -2e- lxp T* f {t,t') - e~ 2tv TZ f {t,t') = Tl{{t - t'). 



(3.19) 



Note that all these commutators are c-number fun ctions and this is essentially a consequence of the commutation 
rules between the noise operato rs giv en by Eqs. (|2.3|) and (2.15). At this point it is possible to compute the correlation 
func tion C (t,t') replacing Eq. ( 3.11 ) in Eq. (|3.9| ): we obtain an integrable differential equation of the same form of 
Eq. (3.10), whose solution is 



C(t,t') = C(0,0) X(t)x(t') + W\a)G(t,t'), 



(3.20) 



where 



Q{t,t') = l 



m\i\[t,t'\ 



dt" x (t - t") x (t' - t") 



/•min[i,£ — r] 

-g sin(0 - ^)0(t' - r) / dt" X (t - t") X (t' -r-t") 

-g sin(0 - ip)S(t - r) / dt" x(t-T- t") X (t' - t") 



(3.21) 



o 2 sin 2 ^- w) /■min[t-T,t'-T] 

+ g l yj 9(t - r)0(t' - r) / dt" x(t-r- t") X (t' - t - t")\. 

V Jo J 

If we now set t 1 = t in this expression, we get the second order moment (Eq. (^^) with ./V = 2); it is however more 
useful to consider the expression of the following "variance" 



a^(t) = 2{ 



(X 2 (t)), a f(X v (t)) 



V (B\a) ) J 



1 + JL g 2 sin 2 ( 

2 2t? V 



vmt-r) dt" x 2 (t")- 

Jo 



(3.22) 
(3.23) 



Let us consider again the physically interesting limit of small delay, 7T <C 1, in which the variance of Eq. (3.23) can 
be well approximated by the first order expansion in 77-, which is given by 



1 



9 

9 ---2, 



2 I r\ 



1 _ e -"/(l-2gsm(6-<p))t 
1 -2gsm(9-ip) 



(3.24) 



g 2 sm 2 (6 - cp) 
2^ 



(l + 7 i#sin(0 - (p))e~' ({1 - 29sin(e ~' fi) '> t + gsm(6 - <p) 



1-e 



- 7 (l-2gsin(e- ¥ ;))i 



1 -2gsin(6-ip) 



JT . 



This expression, as well as Eq. (ETa), is valid for t > r only, since for t < r the feedback is not yet acting and the 
variance assumes its value in absence of feedback, a^ 2 '{t) = 1/2. 

As concerns the higher order variances, it is convenient to consider the following quantities 



/ \ N 

(X t J(t)) l3a = ([X v (t)-x(t)X lp (0)) ) 



(3.25) 



and to proceed as in the previous case, that is, by considering the function (X^ (tyXptt')) p a , and differentiating 
it with respect to t! by keeping t constant. This gives a differential equation which can be formally integrated and 
setting then t = t! it is easy to get the following recursive relation: 



(*£(*)>*, = (N- l)C?(t,t)(l^- 2 (i)), Q . (3.26) 

This relation can be easily solved and it can be expressed in the following way 

{yX v (t) j }p a = {^X tp (t) x{t)) }f> a (3.27) 

for N odd 




(3.28) 

o^{t)] N/2 for N even, 

reproducing the results for the mean and the variance derived above for N — 1 and N = 2 respectively. 



The moments of Eq. (3.28) satisfy the typical relation of a Gaussian process and this provides an independent check 
of the fact that the probability distribution P(x Vl t), being a superposition of Gaussians at t — 0, remains Gaussian 
at all times, as it must be, due to the linearity of the evolution equation. This probability distribution can be written 
in the following form 



^-V.^-im-'t J^~ )• (3 - 29) 



with x(t) and <J^ 2 \t) given respectively by Eq. ( |3.6| ) and Eq. (3.23). If we set t = in these expressions, we obtain 
the exact solution of the ideal case of zero feedback delay , which has been derived in p] . 

It is instructive to apply the general result of Eq. (3.29|) to the case of an initial even Schrodinger-cat state 



\^Cat) T = iVJK) + | - «o)} <8> |0), (3.30) 

in which | ± ao) are two coherent states of the cavity mode and N = (1 + e~ 2 l"°l )~ 1 / 2 is the normalization constant, 
to see the effect of the non-zero delay on the decoherence process. The marginal probability distribution P(x v , t) for 
this initial condition can be written as (see also Ref. J6|) 

P(x v ,t) = N 2 {p 2 + (x v ,t) + p 2 _(x (p ,t) + 2p + (x v , t)p-(x v ,t) coB[n(x v ,t)](ao\ ~ a )" (t) }, (3.31) 

where the first two terms, 

p±M = T^m exp i ^w — | (3 - 32) 

pertain to the two initial coherent state, while the third, containing the functions 

ne>.0 = 3 " tn X" >X( " ( " 3) 

and 

describes the time evolution of the quantum interference between them. In Ref. B it is shown that, in the case of zero 
feedback delay time (r = 0) and perfect homodyne detection r\ = 1, this interference term decays with a decoherence 
time 

tdec{9) ~ 2 7 |a | 2 (l- 5 sin(^)r (3 ' 35) 

which for < gsm(9 — (p) < 2 implies that quantum coherence survives for a longer time with respect to the no- 
feedback case (g = 0) . If we consider the presence of a non-zero delay time, the correction to this feedbac k-indu ced 
decoherence slowing down can be evaluated from the behavior of the fringe visibility function rj(t) of Eq. ( |3.34 ). It 



is however more instructive to see the effects of the feedback delay on the plots of the probability distribution. In 
Fig. g and Fig. || we show the plots of P(x v ,t) (x = x v= o) for the case ao = z*5. What is relevant in Fig. 2 is that 



the probability distribution in the presence of a non-zero delay 7T = 0.01 and efficiency 77 = 1 (Fig. 2(c)) looses its 
interference fringes only slightly faster than the ideal case of zero delay and 77 = 1 (Fig. 2(b)) and that the decoherence 
is still much slower than the no-feedback case (Fig. 2(a)). Moreover in Fig. 2(c) it is once again very evident the "step 
effect" , producing a rapid initial "flattening" of the probability distribution, which is quantitatively the main effect 
of the feedback delay. Fig. 2(d) shows the effect of a non-unit detection efficiency ( 770. 9) which, as it can be easily 
expected, degrades the performance of the homodyne feedback scheme in an appreciable way. 

Fig. 2 and Fig. 3 show that an appreciable decoherence retardation is obtained when the condition 7T < 0.01 is 
satisfied and this means a feedback delay time equal to one half of the decoherence time in absence of feedback, 
tdec(0) — O.O27 -1 . Therefore, the feedback-induced decoherence retardation takes place even in the presence of a 
non-zero delay and Figs. 2 and 3 show that one can even tolerate delay times of the order of the decoherence time 
itself. 

B. Complete solution of the dynamics 

In this section we exactly solve the time evolution of the cavity mode in terms of the symmmetrically ordered 
characteristic function, which is nothing but the expectation value of the cavity mode dispacement operator on the 
initial state of the whole system 

Xw(X,t)=Tr( PT {0)D{\,t)) D{\,t)=exp(\a^t)-\*a(t)). (3.36) 



(3.37) 



Using the commutation rules of this operator with a(t), Eq. (2.16) becomes in this case: 
dD(X,t) = %l\* D{\,t)a{t) - \<J{t)D{\,t)Xdt+ ^\*D{\t)dB{t) 
- /7 XD(X,t)dB j '{t) -SLg 2 (\ e - ie + X*e M ) 2 D(X,t)e(t-T)dt 

+1 ^g(Xe' ie + X*e i6 ) { ^a){t - r) e* dt + dB ^~ r) e ^ } D(X, t) @{t - r) 

+1 ^g(Xe- lt> + X*e lf> ) D{X, t) { ^ja(t - r) e"^ dt + iMtzll e -i<P J Q( t _ r ). 
When we consider the initial condition of Eq. ([3J) , we can simply focus on the matrix element 

(D(X,t)) fia = {(0| ® (P\}D(X,t) {\a) ® |0>}, (3.38) 

which obeys the following evolution equation 

f t (D(X,t)) 0a = %{x*{D(X,t)a(t)) 0a -X{aHt)D(X,t)) 0a ] (3.39) 

-n- g \Xe-^ + X*e^Y{D{\t)) fi Mt-T) 
+q g {Xe- i$ + X*e w ) { e^(at(i _ T ) D{X,t)), a 

+e-^(D(X,t)a(t - T)) (ia }0(t - t). 



To solve Eq. fl3.39| ), we have to deal with terms of the form (D(X,t) a(t'))p a . We first focus on its t' dependence, 
which can be determined in the same way as we have done in the previous Section for (X ¥ ,(i)X I/ ,(i')) /3a , that is, b y 



differentiating with respect to t 1 (D(X,t) a(i')) /3a , by keeping t constant. Using the commutation rules of Eq. (3.15) 
it is possible to derive the following relations, valid for all positive value of t and t': 

[D(X,t),a(0)] = D(X,t)Vx(t) [^(A,i),a f (0)] = D(X,t)V* x (t) 

[D(X,t),dB(t')] =D(X,t)K x (t-t , )dt' [D(X,t),dB*(t')] = D{X,t)K\{t - t')dt' , (3.40) 

[D(X,t),dB f (t')] = D(X,t)K{(t - t')dt' [D(X,t),dBJ(t')] = D(X,t)Tl{*(t-t')dt', 
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where we have defined 



Va(<) 



2 v > v ' sin(6-ip) 



(3.41) 



TZ x (t) = XRi(t)- \*H(t), 



(3.42) 



n{(t) = m f 1 {t)-\*ii f {t), 



(3.43) 



(H(t), Hi(t), TZ f (t), and n{(t) are given by Eqs. d3~16| ), < $TJ\ ), (fl|) and (fig))). Using Eqs. ( |3~40| ), we then obtain 



where 



(£>(A,t)a(t')>^ = (£>(A,t)), Q {ae-**' - ^ e (/3V^ + ae-^)0(t' - r) x( * ' ~ ' 



-it' 



sin(9 — tp) 



+e~* 



f '- T )e(t'-r)W A (t'-r,t)}, 



W x {t> - r,t) = -i^-ge 1 



t'-r 



dtTe* 



'{2^yA x (t",t) 






(3.44) 



(3.45) 



A* (*,*') 



V7 



e iv I dt" X (t-t") ®(t' - t")nX(t' - t") - gsm(8 - tp)Q(f - r) ^ A ^ - + ^ 

Vv 



+^x(tM(t'). 



(3.46) 



Using the fact that (a^(i') D(X, t)) fja = (D(— A, t) a(i'))* /3 and Eq. (3.44), Eq. ( 3.39J ) becomes the simple homogeneous 
differential equation 



-{D(X,t)) 0a =H(t){D(X,t)), a 
whose solution is 

(D{X,t)) 0a = (D(X,0))p a exp f f dt"H{t") 
This result is only apparently simple, since the explicit time dependence of T~L(t) is given by 
H(t) = J (A* a - A/3*)e-i* + Uxe~ i$ + X* e w ) (0* e* + a e - l *)Q(t - r)^? + *' 



$e-** 



sin(# — if) 



-^-g 2 (Xe- ie + X*e l9 ) 2 Q{t - r) + |e(* - t^-^-^ [x*W x (t -T,t) + XW* x {t - r, t) 



+ilg(Xe- te + X*e t6 )Q(t - 2r)e-^*- 2T ) e~ iv Wx(t - 2r, t) - e^W^t - 2r, t) 



(3.47) 



(3.48) 



(3.49) 



Eqs. (3.48) and (3.49) describe the time evolution of the cavity mode starting from the initial condition ( p.2| ), in the 
case of a non-zero feedback delay time. It is however interesting to consider the approximated expression of this result 
at first order in jt since this condition can be easily realized experimentally with usual electro-optical feedback loops 
and good cavities. 
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C. Approximated expression for {D(X,t)) /la in the jt <C 1 limit 



We have two possible equivalent ways to deal w ith the jt <^ 1 limit. The most straightforward one is simply to 
consider this limit in the exact solution Eq. (3.48). However this procedure is not very trasparent from the physical 
point of view because of the complicated form of the function 7i(t). It i s inst ead more instructive to perform the same 
limit from the beginning on the evolution equation for (D(X, t)) ?a , Eq. (3.39), and then integrate it. We shall consider 
this second approach also because it can be adopted not only in the problem considered in this paper (linear choice 
(3.1) fo r the feedback operator F(i)) but also for more general forms of the operator F(t). Let us go back therefore 
to Eq. ( 3.39 ), where the difficult terms to handle are those containing (D(X, t) a(t — T)) 0a and its complex conjugate, 
which can be rewritten as 



where 



(D(X,t)a(t - r)) 0a = (D(X,t) a(t)) 0a - (D(X,t) Aa(t - T )) Pa , 



Aa(t -r) = a(t) - a(t - r). 



(3.50) 



(3.51) 



In the limit jt -C 1, we can use Eq. (2.16) to approximate Eq. ( |3.5l| ) at the first order in r so that we can write (we 
also consider t > 2r) 



(D(X, t) a(t - r)) ffa ~ (D(X, t) |o(t) + %a(t-T)T + i 91 e l6 X v (t - 2r) r}), : 

+ i9^ke i{6+tp) (D(X,t) AB}(t-2r)) flal 
where ABf(t — 2r) is the following Ito increment 

AB f (t - 2r) =B f (t-r)- B f (t - 2r). 



The last term in the right side of Eq. ( |3.52| ) can be simplified using the identity 

D(X, t) = D(X, t - 2r) + { D(X, t) - D(X, t - 2r) } 
and approximating the term in the curly brackets again at first order in r. Finally one gets 

(D(A, t) a{t - r)) 0a ~ (l + ^(1 + ige l{6 -^) ) (D(X, t) a{t)), a 



(3.52) 



(3.53) 



(3.54) 



(3.55) 



+ i 71 e i(o+ v ) {a i {t)D{X!t))0o 



^-gi(Xe~ w + X*e l9 y e (D(X,t)) lia . 



Replacing this approximated expression in Eq. ( |3.39| ) together with the corresponding one for (a^(t — r) D(X,t)) f3a , 
one obtains a simple integrable equation (valid for t > 2r) having the following solution: 



(D(X, t)) 0a = (fi\a) exp |-^(i)|A| 2 + ^(^A 2 + B 2 {t) (A*) 2 + C(t)X + V(t)X* } 



where 



A{t)= l -+ 9 -{Mt)- 1 rT 1 (t)} 

Bi(*) = (fi 2 (t))* = -|^{^ (i)-7r^iW}e- 219 , 



in which 



M*) 



l-e 



-(l-2gsine)7* 



1 - 2g sin6» 



T^t) = (1 + 7 ;.gsin#)e-( 1 - 2 9 sin ^ f + (g am0)T o (t), 



(3.56) 

(3.57) 
(3.58) 

(3.59) 
(3.60) 
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and 



C{t) 



-io 



2sin6» 



( oe -« + /3VVi* + ^-5(0 + /3*)e-( 1 - 2 " fa0 ^ t 



2 sin6» 



+ |i 5 e- ie (a +/T)((1 -2 gS m0)|t-l) e -( 1 - 2 ^ ine )* t | 7 T, 



(3.61) 



Z>(t) = - 



— (ae~ ie + pe^e-l* + ^-5 (a + ^-(i-a***)^ 



2sin# 



2 sin6» 



-,ge 



i6 {a +/?*)((! 



2#sin6»)-£ 



le 



-(l-2gsin9)2.t 



7T. 



(3.62) 



(We have chosen the phases so t hat 7 = 0). As expected, setting t = one obtains the same results of Ref. ||. 
It is important to note that Eq. ( 3.56 ) has been obtained integrating from t = 2r, i .e., using (D(A, 2r)) /Sa as initial 
condition. This initi al con dition could be derived from the exact solution Eq. ( 3.4S| ), but it could also be obtained 
by noting that Eq. ( 3.39 ) takes a very simple form for < t < 2r. In fact the terms (D(X,i)a(t — r)) 0a and 
(a* (t — t) D(X, t))p a , for r < t < 2r, can be written as 



-i(t-T) 



(D{\, t) a(t - r)) 0a = a(D(X, t)) 0a e'^ 



{a\t - t) D(X, t)) 0a = 13* (D(X, t)), a e-i^ T \ 



(3.63) 



depending on the fact that for r < t < 2r the operators a(t — t) and a) (t — r) are not affected by the feedback loop. 
Moreover, for < t < r these ter ms d o not contribute to the evolution because of the presence of the step function 
in the equation. In this way Eq. ( 3. 39] ) can be simply integrated for < t < 2r too. 

In Ref. |6| , the decoherence inhibition capabilities of the homodyne-mediated feedback scheme have been described 
by looking at the so-called coherence function C(t) = (—ao\p(t)\a>o) ■ However, an equivalent description of the 
decoherence of the Schrodinger- cat s tate (|3.30[) is provided by the time evolution of the expectation value of the 



operator D(2ao,t) on the state ( 3.30J ). In fact, for |ao| 2 3> 1 and jt <C 1, this quantity has the same time behavior 



of the coherent function C(t) of Ref. J6), because its off-diagonal contributions are greater then the diagonal ones by 
the exponential factor e 2 '" ' . In the small delay time limit 71c 1, this equivalent coherence function takes a simple 
form, given by 



(£(2a ,;))^iexp{-2|ao| 2 [2 + i^fl^ (t) - 7 t ^(t)) 



-2e-i^- 29sine ^ + (gsm0)[2 - jt(l - 2gsm9)] 1 Te-K 1 - 



2gsin8)t 



Y 



(3.64) 



where we have defined phases so that ip — and considered Re{ao} = 0. I n the zero delay time case (r = 0), 
whe n jt < C 1, this function is an exponential with characteristic time equal to ( 3.35J ) (see Ref. raj). We have plotted 
Eq. ( 3.64 ) in Fig. ^ again for ao = i5 in the ideal case rj = 1 for different values of the feedback delay and we have 
compared it with the no- feedback case (dot-dashed line) and with the ideal case r = and r/ = 1 (dashed line). As we 
have seen above, the decoherence slowing down is significant up to jt = 0.01, i.e. r = Q.btdec and of course increases 
as long as the feedback delay decreases. In Fig. 5 the effect of an imperfect homodyne detection is considered and 
it is shown how, as expected, the decay of the coherence function D(2ao,t) becomes faster and faster as long as the 
efficiency 77 decreases. 



IV. CONCLUSIONS 

In this paper the dynamics of a cavity mode subject to a feedback loop using a fraction of the output homodyne 
photocurrent to control the transmittivity of a mirror has been completely solved in the general non-Markovian case 
in which the feedback delay time is not negligible. To solve this problem we have generalized the derivation of the 
feedback quantum Langevin equation based on the input-output theory of Ref. 0] to the non-unit detection efficiency 
case. We have seen that when rj < 1, the fed-back photocurrent involves a new Ito noise term, coinciding with the 
usual input noise only in the 77 = 1 limit. 
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We have seen that the main effect of the delay is the "step effect", that is, the fact that the feedback loop starts 
acting only after the time r has elapsed. Apart from this, the dynamical behavior in the presence of a feedback delay 
does not differ very much from the predictions of the Markovian treatment based on the zero-delay limit, as long as 
7T -C 1. As a consequence, the significant decoherence slowing down demonstrated in Ref. |6j for the zero delay limit 
holds even in the presence of a non-zero delay r, and one gets a good decoherence inhibition for values of r up to one 
half of the decoherence time. 
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APPENDIX A: 

Let us consider the evolution equation for the cavity mode state vector \4>w{t)) for the case of feedback with zero 
delay time, given in Q6Q : 

d\i> w (t)) = {Adt + Bdw(t)}\ip w (t)), (Al) 

in which dw(t) is a real-valued Wiener increment deriving from the input noise, and 

A = -^a-^F 2 -i 1 Fae- llp , (A2) 

B = V7(ae"^ - iF), (A3) 



are cavity mode operator s in the Schrodinger picture, with F = g(ae +a'e )/2 as in Eq. (3.1): this equation 



can be derived from Eq. ( 2.1£ ) by means of the method developed in [[L9| . Now we adopt the stochastic integration 



method derived in pll to cxplicitely solve Eq. (Al). The explicit analytical solution is given by 

\ip w (t))=E w (t)\iP w (0)), (A4) 

where the evolution operator is 

B 2 
E w (t) = exp [iA(t)] exp [(A - — ) t] exp [x-(t)a + x+(t)^] , (A5) 

with A(t) and X±(t) complex functionals of the Wiener process w(t); defining in particular the functions 

M) = ^{~ igil+e-^e^e^ + [- ige- ie + 2e"^ + zge^'^e'^} (A6) 

/*(*) = ^{- SO-' e-^yOel* +%[- ige~* e + 2e~^ + ige^' 2 ^]e'^} (A7) 



we have that 



Xj(t) - / fj(t')dw(t') for j = l,2 (A8) 

Jo 

Ai/)= J h(t') X 2{t')dw{t')- J' h{t') X i(t')dw(t') (A9) 

X±(*) = -^(xi(*)±*Xa(*)). (A10) 



If we consider the case of an initial coherent state |ao), Eq. flAq ) gives 
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|^(*)>=fi«,(*)|[x+(t) + oo]e-**) (All) 

that is, the state remains coherent with amplitude [%+(£) + a o]e~^S and £ w {t) is a complex weight given by 

£ w {t) = exp{iA(*) + !(x+(*)* + X-(*))(X+(*) + 2 «o) + i Im[x+(*K] (A12) 

+ jf TfeW-ri - | [| X+ (t) + a \ 2 + ( x+ (t) + a ) 2 e- 2 ^] (1 - e"^)}. 



First of all we note that replacing g = in Eq. (All) we correctly obtain the same results of Ref. Jj22| for the case 
of no feedback loop. The state remains coherent in agreement with the "no-go" theorem of Ref. 07 according to 
which homodyne- mediated feedback is not able to increase the nonclassicality of the output light. Using linearity, 
from Eq.( All ), one easily gets the evolution of the Schrodinger-cat state 



\ipCat)=N\\a ) + |-ao}} 

and it is also possible to evaluate the coherent function defined in H, 

-a \^ w (t))(ip w (t)\a ) 



C w (t) 



ftM*)|iM*)> 



(A13) 



(A14) 



By considering the limits |«o| ^> 1 and yt <C 1 and setting set (p = and Re{«o} = as in ||, one gets the following 
expression 



C w {t) = - exp{ - -yg 2 w 2 (t) - 2*V7|ao|(l - 9 sin6)w(t)], 



(A15) 



reproducing very well the numerically obtained stochastic trajectories of Ref. q|. Eq. ( A15 ) clearly shows the results 
of Ref. M, i.e. that the modulus of C w (t) is not (9-dependent and that the large fluctuations in the absence of feedback 
and those in the presence of feedback but with phase 9 ^ ir/2 are essentially phase fluctuations. 



APPENDIX B: 



Let us consider the equation (|3.4|) for the function (X v (t)) AB , 

d 



dt 



Z(0 = -Z(0 + 2kZ(C-y)e(C-y), 



(Bl) 



where for sake of simplicity we have introduced the following notation: 



f — 21 
** 2 ' 



'YT 



z(0 = (x v ( 2 4)h 



k = gsm(9 — ip). 



(B2) 



We are interested to study Eq. (Bl) for £ > 0, with initial condition Z(0) = Zq. First of all we observe that Z(£) 
is a continuous function of class C°°(I) on every interval I = 1 ny, (n + l)y [ , with n G N. Moreover it is possible 
to verify that its p th derivative is not cont inuous at £ = py. Apart from this irregular behaviour due to the presence 
of 0(£ — y), it is easy to see from Eq. (Bl) that the generic solution Z(£) is bounded by a locally integrable function 
and therefore it can be Laplace-transformed. Denoting with Z{s) the Laplace transform of Z(£), from Eq. ( p31| ) we 
have 



Z(s) 



Z Q 



(* + l)(l 



2ke~sy > 



(B3) 



It is always possible to choose the integration path for the antitransformation so that 
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2ke~ sy 
s+1 



<1, 



(B4) 



and using the geometrical series, we can write (B3) as 



+00 



Z(s) = Z J2 



(2k) n e 



sny 



(s + 1)"+! 



(B5) 



which can be easily antitransformed. Using relations Eq. (B2) we then obtain the solution (3.5) of Section 2.2 
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FIG. 1. Normalized mean value of the measured quadrature \(t) (see Eq. (3.5)) for gsin(9 — ip) = (without feedback, 
dot-dashed line), gsin(8 — ip) — 0.45 and some values of the delay time: from bottom to top: 7T = 5, 7T = 2.5, 7T = 1, 
"/t = 0.5 (full lines) and 'yr — (dashed line). 
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FIG. 2. Time evolution of the marginal probability distribution P(x,t) (x = x v ) for the initial state ( 3.3C ) with Qo = »5. 
(a) refers to the no feedback case (gsin(9 — ip) — 0); (b) refers to the ideal limit of zero feedback delay time and r\ — 1, with 
gsin(8 — ip) — 1; (c) refers to the case with a non-zero feedback delay time (-yr — 0.01) r) — 1 and with gsin(6 — ip) = 1; (d) 
refers to the case 7T = 0.01, ?7 = 0.9 and <jsin(# — ^>) = 1. 
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FIG. 3. Comparison at -yt = 0.1 of the marginal probability distribution P(x,t) between the cases of no feedback (dot-dashed 
line), feedback with zero delay time (gsin(6 — tp) = 1, r\ — 1, dashed line), feedback with non-zero delay time yr = 0.01 and 
7] — 1 (gsin(9 — (f) — 1, full line), and feedback with non-zero delay time and imperfect detection 7T = 0.01 and r] = 0.9 
(gsin(9 — ip) — 1, dotted line) 
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FIG. 4. Time evolution of 2{D(2ao,t)} for ao = *5 and ^ = for gsinO — (dot-dashed line), gsin9 — 1 with 7T = (dashed 
line) and gsin# = 1 with, from bottom to top, -yr = 0.02, yr = 0.01, 7T = 0.001 (full lines). 
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FIG. 5. Time evolution of 2{D(2ao, £)) for ao = i5 and ip = in the presence of feedback with non-zero delay time (jt = 0.01 
and gsin9 — 1) and for different values of the homodyne detection efficiency r\\ from bottom to top, r\ = 0.75, r\ = 0.9, r\ = 0.95, 
7] — 1 (full lines). The dashed line refers to the no-feedback case. 
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